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The impact of cluster mergers on arc statistics 
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ABSTRACT 

We study the impact of merger events on the strong lensing properties of galaxy clusters. Pre- 
vious lensing simulations were not able to resolve dynamical time scales of cluster lenses, 
which arise on time scales which are of order a Gyr. In this case study, we first describe qual- 
itatively with an analytic model how some of the lensing properties of clusters are expected 
to change during merging events. We then analyse a numerically simulated lens model for 
the variation in its efficiency for producing both tangential and radial arcs while a massive 
substructure falls onto the main cluster body. We find that: (1) during the merger, the shape of 
the critical lines and caustics changes substantially; (2) the lensing cross sections for long and 
thin arcs can grow by one order of magnitude and reach their maxima when the extent of the 
critical curves is largest; (3) the cross section for radial arcs also grows, but the cluster can effi- 
ciently produce this kind of arcs only while the merging substructure crosses the main cluster 
centre; (4) while the arc cross sections pass through their maxima as the merger proceeds, the 
cluster's X-ray emission increases by a factor of ~ 5. Thus, we conclude that accounting for 
these dynamical processes is very important for arc statistics studies. In particular, they may 
provide a possible explanation for the arc statistics problem. 



1 INTRODUCTION 



The abundance of strong gravitational lensing events produced by 
galaxy clusters is determined by several factors. Since light deflec- 
tion depends on the distances between observer, lens and source, 
gravitational lensing effects depend on the geometrical properties 
of the universe. On the other hand, gravitational arcs are rare events 
caused by a highly nonlinear effect in cluster cores. They are thus 
not only sensitive to the number density of galaxy clusters, but also 
to their individual internal structure. Since these factors depend on 
cosmology and in particular on the present value of the matter den- 
sity parameter flrjM an d the contribution from the cosmological 
constant floA, arc statistics establishes a highly sensitive link be- 
tween cosmology and our understanding of cluster formation. 

Using the ray-tracing technique for studying the lensing 
properties of galaxy clusters obtained from N-body simulations, 
Bartelmann et alJ 119981) showed that the number of giant arcs, 
commonly defined as arcs with length-to- width ratio larger than 10 
and fi-magnitude smaller than 21.5 fWul 19931) . which is expected 
to be detectable on the whole sky, differs by orders-of-magnitudes 
between high- and low-density universes, strongly depending even 
on the cosmological constant. In particular, they estimated that the 
number of such arcs in a ACDM cosmological model (Qqm = 0.3, 
^0A = 0-7) is of the order of some hundreds on the whole sky, 
while roughly ten times more arcs are expected in an OCDM cos- 
mological model (£2om = 0-3, Qqa = 0). 

Although still based on limited samples of galaxy clusters, ob- 
servations indicate that the occurrence of strong lensing events on 



the sky is rather high jLuppino et allll999l : IZaritskv & Gonzalez! 
l2003tlGladders et all2003T) . For example, searching for giant arcs 
in a sample of X-ray luminous clusters se lected from the Ein- 
stein Observatory Extended Medium Survev. lLuppino et alJll999T) 
found that their frequency is about 0.2 — 0.4 arcs per massive clus- 
ter. Despite the obvious uncertainties in the observations, the only 
cosmological model for which the number of giant arcs expected 
from numerical simulations of gravitational lensing comes near 
the observed number is the OCDM model. In particular the ob- 
served incidence of strongly lensed galaxies exceeds the predic- 
tions of a ACDM model by about a factor of ten. On the other hand, 
based on the observation s of type la supernovae IPerlmutter et alJ 
Il999t lTonrvetalJl2003l) and the recent accurate measurements 
of the temperature fluctuations of the Cosmic Microwave Back- 



ground obtained with balloon experiments (e.g. Istompor et 



2001; laffe et al. 2001; Abroe et al. 2002; Benoit et al. 2003) or by 



the WMAP mission (e.g. lBennett et alJl2003HSpergel et alJ 



al. 



2003), 



the ACDM model has become the favourite cosmogony. This is 
known as the arc statistics problem: the mismatch between the 
observed number of arcs and the number expected in the ACDM 
model preferred by virtually all other cosmological experiments 
hints at a lack of understanding of cluster formation. 

Several extensions and improvements of the numerical sim- 
ulations failed in finding a solution to this problem in the lens- 
ing simulations. Meneghetti et al. ( 2000) studied the influence of 
individual cluster galaxies on the ability of clusters to form large 
gravitational arcs, finding that their ef fect is statistica lly negligible. 
ICooravl 1 19991) and Kaufm ann & Straumannl fcOOOl) . using spher- 
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ical analytic models and the Press-Schechter formalism for mod- 
elling the lens cluster po pulation, predicted a weaker dependence of 
arc statistics on £2oa, but Meneg hetti et alJ l2003bl) . comparing nu- 
merical models of galaxy clusters and their analytical fits, showed 
that analytic models are inadequat e for qu antitative studies of arc 
statistics. Finallv. lMeneghetti et al.li2()0 3a1 found that the presence 
of central cD galaxies can increase the cluster efficiency for pro- 
ducing giant arcs by not more than a factor of about two. 

As alternat ive solution of the arc statistics problem, 

iBartelmann et alJ 120031) recently investigated arc properties in 
cosmological models with more general forms of dark en- 
ergy than a cosmological constant. Several studies showed that 
haloes should be more concentrated in these models than in 
cosmological-constant m odels with the same dark energy d ensity 
today iBartelmann et alJl20o3: iMainini et ail 120031 ; iKlypin et all 
120031 iDolag et alj|2003h . allowing them to be more efficient for 
strong len sing. U sing simple models with constant quintessence 
parameter, Bartelmann et al] 120031) found that the relative change 
of the halo concentration is not sufficient to produce an incre- 
ment of one order of magnitude in the expected number of giant 
arcs. Nonetheless, other more elaborate dark-energy models need 
to be evaluated numerically and will be discussed in a future paper 
(Meneghetti et al. in preparation). 

In this paper we investigate another possible effect which 
could not be properly considered in the previously mentioned nu- 
merical simulations of gravitational lensing by galaxy clusters. In 
those works, the lensing cross sections for giant arcs of each numer- 
ical model were evaluated at different redshifts, with a typical time 
separation between two consecutive simulation outputs of approx- 
imately At ~ 1 Gyr. Therefore all the dynamical processes arising 
in the lenses on time scales smaller than At were not resolved. 

N-body simulations show that dark matter haloes of differ- 
ent m asses continuously fall onto rich clusters of galaxies iTormenl 
1997). The typical time scale for such events is ~ 1 — 2 Gyr, which 
therefore might be too short for having been properly accounted for 
in the previous lensing simulations. 

However, the effects of mergers on the lensing properties of 
galax y clusters may pote ntiall y be very important. As d iscussed by 
Bartelman n et alJ 1 19951) and Menegh etti et alJ <2003bl) . substruc- 
tures play a very important role for determining the cluster effi- 
ciency for lensing. Indeed, analytic models, where substructures 
and asymmetries in the lensing mass distributions are not properly 
taken into account, systematically underestimate the lensing cross 
sections of the numerical models. The main reason is that mass 
concentrations around and within clusters enhance the shear field, 
increasing the length of the critical curves, and consequently the 
probability of forming long arcs becomes higher. 

Given the strong impact of substructures on the lensing prop- 
erties of galaxy clusters, it is reasonable to expect that during the 
passage of a massive mass concentration through or near the cluster 
centre, the lensing efficiency might sensitively fluctuate. While the 
substructure is approaching the main cluster clump, the intensity of 
the shear field and, consequently, the shape of the critical curves 
may substantially change. Moreover, while the infalling dark mat- 
ter halo gets closer to the cluster centre, the projected surface den- 
sity increases, making the cluster much more efficient for strong 
lensing. 

This paper describes a case study on how much the lensing 
cross sections change during the infall of a massive dark matter 
halo on the main cluster progenitor. For doing so, we investigate the 
lensing properties of a numerically simulated galaxy cluster during 
a redshift interval when a major merging event occurs. The gen- 



eral aim is to understand if mergers can enhance the cluster lens- 
ing efficiency sufficiently to provide a solution to the arc statistics 
problem. 

The plan of the paper is as follows. In Section 2 we use a 
simple analytic model based on the NFW density profile for com- 
puting the strong lensing properties produced by cluster mergers. 
In Section 3 we describe our lensing simulations utilizing a nu- 
merical model obtained from a high-resolution N-body simulation; 
a method for evaluating the dynamical state of the cluster is also 
introduced. In Section 4 we present our results for critical lines 
and caustics, and we estimate the expected number of tangential 
and radial arcs. Section 5 is devoted to a discussion of the obser- 
vational implications of our results, and we summarize the results 
and present our conclusions in Sect. 6. 



2 ANALYTIC MODELS 

We begin by investigating with the help of analytic models how the 
lensing properties of a cluster lens are expected to change while a 
massive substructure passes through its centre. We will model both 
the main cluster clump and the infalling substructure as spherical 
bodies with the density profile found in numerical simulations by 
lNavarroetai]ll996l hereafter NFW). 



2.1 NFW model 

The NFW density profile is given by 

pW 



(r/r,)(l + r/r,)2 



(1) 



where p s and r s are the characteristic density and the distance scale, 
respectively (see Nav arro et all 19971) . The logarithmic slope of this 
density profile changes from —1 at the centre to —3 at large radii. 

We briefly summarize here the main features of NFW haloes. 
First, p s and r s are not independent. Second, the ratio between the 
radius tiqq, within which the mean halo density is 200 times the 
critical density, and r s defines the halo concentration, c = r2oo / '"s • 
Numerical simulations show that this parameter systematically 
changes with the halo virial mass, which is thus the only free pa- 
rameter. Several algorithms exist for computing the concentration 
parameter from the halo v irial mass iNavarro et alJl99"7l : lEke et alJ 
2001; Bullock et al. 2001). In the following analysis we adopt the 
method proposed by Nava rro et al 1 <1997t) . Third, the concentra- 
tion also depends on the cosmological model, implying that the 
lensing properties of haloes with identical mass are different in dif- 
ferent cosmological models if they are modelled as NFW spheres. 



discussed in several papers (Bartelmann 1996 


Wright & Brainerd 


2000; Wvithe et al. 2001; 


Li&Ostriker 2002; 


Perrottaet al. 2002; 


Bartelmann et al. 2002; Meneghetti et al 







Due to axial symmetry, lensing by an NFW sphere reduces to a 
one-dimensional problem. We define the optical axis as the straight 
line passing through the observer and the centre of the lens, and 
introduce the physical distances perpendicular to the optical axis 
on the lens and source planes, % and T|, respectively. We introduce 
dimension-less coordinates x = ^/^o an d y = T l/ T lo by means of 
the coordinate scales i;o = r a and "Ho = {D s / D\)v\q, where D s and 
D\ are the angular diameter distances to the source and lens planes, 
respectively. 

Using this notation, the density profile implies the surface 
mass density I B artelmanrl 1 99d) 
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Figure 1. Illustration of the analytic model described in Sect. 12. 2| Maps of the tangential eigenvalue X, for all configurations of the main cluster body relative 
to the infalling substructure. The critical lines are over-plotted in each map. The (comoving) size of each panel is ~ 1.2/i Mpc. Different panels refer to 
different distances d (in units of the scale radius of the larger halo) between the two haloes, as indicated by the top label. 
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The lensing potential is given by 

>r(x) = 4K s s(*), 



where 



8W 



2*2 



2 arctan 2 
-2arctanh 2 



(x>l) 

(x<l) 

(x=l) 



(2) 



(3) 



(4) 



(5) 



and K s = p s r s I.~ l . here £ cr = [c 2 /(4jtG)] [D s /(DiD ls )] is the crit- 
ical surface mass density for strong lensing and D\ s is the angular 
diameter distance between lens and source. This implies the deflec- 
tion angle 



a(x) - 
with 

h(x) -- 



dx 



4k, 



h(x) , 



ln- + < 



. arctan 1 



= arctanh 1 



(x>l) 

(x<l) 
(x=l) 



(6) 



(7) 



2.2 Expectations 

In this analytic computation, the main cluster body has a mass of 
M ma i n = 10 15 h~ l MQ. It merges with a massive substructure of 
mass M su b = M ma i n /4 = 2.5 x lO 14 /z _1 M , which is a fairly typical 
ratio for merging halo masses. In this simple test, we keep the lens 
plane at the fixed redshift z\ = 0.3, while the source plane is placed 
at redshift z s = 1. We assume a background ACDM cosmological 
model, with SIqm = 0.3, f2(>A = 0.7, and Hubble constant (in units 
of 100 km/s/Mpc) h = 0.7. 

We simulate the infall of the substructure onto the main clump, 
starting from the initial configuration when the smaller halo is 
placed at distance r s from the larger one, where r s is the scale ra- 
dius of the most massive clump. For a halo of mass 1O 15 /i~'M0 
in a ACDM model, r s corresponds to ~ 310 /z~'kpc. Then, we re- 
duce the distance between the two haloes by moving the smaller 
towards the larger one by 0.2r 5 per time step, until the two haloes 
are exactly aligned. 

For any new configuration, we trace a bundle of 512 x 512 
light rays through a regular grid on the lens plane which covers a 
region whose side length is 4r s . This is large enough to encompass 
the critical lines of both the clumps. 

Using we compute the contributions a maln and 6t sub of 
both the main clump and the substructure to the deflection angle of 
each ray. Being linear in mass, the total deflection angle of a ray 
passing through the mass distribution is the sum of the contribu- 
tions from each mass element of the deflector. Thus, the deflection 
angle of each ray parametrized by its grid indices on the lens 
plane is given by 

(I, ■ (i;f . (8) 

The convergence and the shear maps of the lens system are 
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Figure 2. Profiles of convergence K (left panel) and shear y (right panel) along the axis connecting the centres of the main cluster body and the infalling 
substructure illustrated in Fig.Q Different curves refer to different distances d (in units of the scale radius of the larger halo) between the two haloes. Notice 
the increased shear level at fixed x during the merger. 
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Figure 3. Caustics for all configurations of the main cluster body relative to the infalling substructure shown in Fig.Q The size of each panel is ~ 200 h 'kpc 
on a side. Different panels refer to different distances d (in units of the scale radius of the larger halo) between the two haloes. 



reconstructed from the deflection angles. The true position of the 
source on the source plane, y, and its observed position on the lens 
plane, x, are related by the lens equation 



y = x— a(x) 



(9) 



The local properties of the lens mapping are described by the Jaco- 
bian matrix of {5), 



<tyh 



->hk ' 



(10) 



where 8m is the Rronecker symbol. The convergence K and the 
shear components yi and j2 are found from Am through the stan- 



dard relations 



k(x) 



[A n {x)+A 22 {x)\ 



Jl(x) = -~[Ai 2 (x)+A 2 i(x)] 



(11) 
(12) 
(13) 



The Jacobian matrix is symmetric and can be diagonalised. Its two 
eigenvalues are 



Xi — 1 — k — y and 7^- — 1 — k + y . 



(14) 



where y= \/Y? Radial and tangential critical lines are located 
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Figure 4. Left panel: caustic area and length (solid and dashed lines, respectively) as functions of the distance between the merging haloes (in units of the 
scale radius of the more massive halo). The lens redshift is z\ = 0.3. Thick and thin lines refer to source redshifts z s = 1 and z s = 2, respectively. Curves are 
normalized to the maximum value of area and length for sources at z s = 2. Right panel: as in the left panel, but for lenses at redshift Z\ — 0.8 



where the conditions 

Xt = and Xr = (15) 

are satisfied, respectively. The corresponding caustics in the source 
plane, close to which sources are imaged as large arcs, are obtained 
by applying the lens equation to the critical curves. 

We are particularly interested in understanding by how much 
the cluster efficiency for producing tangential arcs changes during 
the merger phase. The size of the lensing cross sections for the 
formation of tangential arcs is strictly connected to the extent of 
the tangential critical curves and caustics: the longer the critical 
curves and the corresponding caustics are, the larger are the lens- 
ing cross sections. Therefore, we first investigate the shape of the 
tangential critical lines and caustics. In Fig.Q we show the maps 
of the tangential eigenvalue A. t for all configurations of the main 
clump relative to the infalling substructure. In each map, we over- 
plot the tangential critical lines. As the two haloes approach each 
other, their critical lines widen and stretch along the approaching 
direction. When the distance is ~ 0.6r s , the critical lines touch and 
merge to form a single critical line. While the substructure overlaps 
with the main cluster body, it shrinks in the direction along which 
the merger proceeds, while it widens perpendicular to it. 

We can explain the evolution of the critical lines as follows. 
As the two haloes approach each other, the shear in the region be- 
tween them grows. This stretches the critical lines in the direction 
along which the two haloes are approaching. When the critical lines 
merge, the shear continues growing in the region inside the unique 
critical line. At this point the critical curve stops growing in a priv- 
ileged direction. It shrinks only in the direction along which the 
haloes merge because of the decreasing distance between their cen- 
tres. On the other hand, it expands almost isotropically by the effect 
of the increasing convergence. We show in Fig.|2|how the conver- 
gence and the shear evolve along the axis connecting the halo cen- 
tres, where shrinking and stretching effects are most evident. The 
main clump is kept at x = 0. Different curves show the results for 
different distances d between the two haloes, in units of the scale 
radius of the largest clump. 

The caustics are shown in Fig. [5] Their evolution reflects that 
of the critical lines. As the shear grows between the two clumps, 
they are stretched and develop cusps. Then, after they merge, the 



single caustic shrinks and reduces to a point when the two haloes 
are exactly aligned. 

In the case of merging haloes, where very complicated struc- 
tures develop in the caustics, it is difficult to say whether the cross 
sections for long and thin arcs scale as the area enclosed by the 
caustics or as the caustic length. In Fig.[4]we show the evolution of 
both the area (solid lines) and the length (dashed lines) of the caus- 
tics while the merger proceeds. Thick lines in the left panel refer to 
the configuration of lenses and sources described earlier. For com- 
pleteness, in the same panel we also show the corresponding curves 
for sources at redshift z s = 2 (thin lines). In the right panel, the cor- 
responding curves for lenses at z\ = 0.8 are shown. As expected, the 
size of the caustics is larger for sources more distant from the lens 
because of its increased convergence. All curves in each panel are 
normalized to the maximum value of length and enclosed area for 
sources at redshift z s = 2. Of course the area enclosed by the caus- 
tics grows more steeply compared to the length of the caustics. It 
is interesting that caustics reach their maximum length before they 
reach their maximum area. This is due to the stretching of the caus- 
tics of the two individual merging haloes towards each other: the 
cusps become very peaked and pronounced until they touch, with 
large change of length. The maximum value of the area enclosed 
is reached only shortly thereafter, when the caustics have already 
started to enlarge isotropically. 



3 NUMERICAL METHODS 
3.1 Numerical model 

We now turn to a numerical cluster model in order to quantify how 
the strong lensing efficiency changes during a merger of a realis- 
tic cluster mass distributions embedded into the tidal field of the 
surrounding large-scale structure. 

The numerical model we investigate here is part of a set 
of 17 objects obtained using the technique of re-simulating at 
higher resolution a patch of a pre-existing cosmologic al sim ula- 
tion. The re-simulation method is described in (Tonne nl ll997t) . A 
detailed discussion of the dynamical properties of the whole sam- 
ple of these simulated cl usters is presented elsewhere iTormen et alj 
l2003tlRasia et all2003l) . 
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Figure 5. Surface density maps of the numerically simulated cluster at several redshifts between z = 0.230 and z = 0.180, during the merger. The scale of each 
panel is 3 ft^'Mpc. 



The parent N-body simulation, with 512 particles in a box of 
side 4797i~' Mpc, h as been produced b y N. Yoshida for the Virgo 
Consortium (see also Jenkins et al. 1998). It assumes a flat universe 
with O-om — 0-3, fi()A = 0-7 and h = 0.7. The initial conditions cor- 
respond to a cold dark matter power spectrum normalised to have 
rjg = 0.9 today. The particle mass is 6.86 x 1O 1o /i _1 M0, which 
allows to resolve a cluster-sized halo by several thousand parti- 
cles; the gravitational softening is 1 5/z 1 kpc. From the final output 
of this simulation we randomly extracted some spherical regions 
containing a cluster-sized dark matter halo. Each of these regions 
was re-sampled to build new initial conditions for a higher num- 
ber of particles - on average 10 6 dark matter particles and the same 
number of gas particles. The mass of the dark matter particles is 
~5x 10 9 A^'Mq, while the mass of the gas particles is ~ 10 times 
smaller. The softening length for dark matter and gas is 3.6 A~'kpc 
and 7.1 /!~'kpc, respectively. These initial co nditions were evolved 
using the publicly available code GADGET iSprinael et aTll200lh 
from a starting redshift z; n = 35 — 50 to redshift zero. 

In our sample of 17 objects we selected a simulated clus- 
ter undergoing a major merger at z ~ 0.2. At redshift z ~ 0.25, 
when their viral regions merge, the main cluster clump and the in- 
falling substructure have virial masses of ~ 7 x m u hr 1 M (3 and 
~3x IO^/j^'M©, respectively. In order to have a very good time 
resolution for resolving the complete merger in detail, we decided 
to re-simulate the cluster between z = 0.25 and z = 0.15, obtain- 
ing 101 equidistant outputs which we use for our following lensing 
analysis. 



ing) side length of the displayed region is 3 A -1 Mpc. We quantify 
the dynamical state of the cluster using the multipole expansion 
techniq ue of the surface density field discussed by Meneahetti et al. 
(2003b). Briefly, we place the origin of our reference frame on 
the centre of the main cluster clump. Starting from the particle 
positions in the numerical simulations, we compute the surface 
density £ at discrete radii r„ and position angles fa taken from 
[0, 1.5] h Mpc and [0,2jt], respectively. For any r„, each discrete 
sample of data T,(r„,fa) is expanded into a Fourier series in the 
position angle, 



E(r n ,^) = ^S,(r„)e- i ^, 
/=o 

where the coefficients S/(r„) are given by 



S,(r„) = £ Z(r n ,fa)d 
k=0 



1% 



(16) 



(17) 



and can be computed using fast-Fourier techniques. 

We define the power spectrum P n (l) of the multip ole ex- 
pansion I as P n {l) = |S/(r„)| 2 . As discussed by Meneahetti et al. 
(2003b), the amount of substructure and the degree of asymmetry 
in the mass distributions of the numerically simulated cluster at any 
distance r n from the main clump can be quantified by defining an 
integrated power Pi nt (r„) as the sum of the power spectral densities 
over all multipoles, from which we subtract the monopole and the 
quadrupole in order to remove the axially symmetric and elliptical 
contributions, 



3.2 Dynamical state of the cluster 

We show in Fig. [5] some projected density maps of the simulated 
cluster at several redshifts during the merger. As is easily seen, a 
large substructure crosses the cluster centre at z ~ 0.2. The (comov- 



%t(r„) = £ />,,(/) -/MO) -/>„(2). 
/=() 



(18) 



This quantity measures the deviation from an elliptical distribution 
of the surface mass density at a given distance r„ from the cluster 
centre. 
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Figure 6. Total integrated power P mi as function of redshift for five different distances r (in units of h Mpc) from the centre of the main cluster clump, as 
indicated in each panel. 



The results are shown in Fig. [6] the total integrated power is 
shown as a function of redshift for five different distances r from 
the centre of the main cluster body. The peaks in each plot arise 
when the infalling substructure enters or exits a circle of radius r. 
Therefore, through the location of the peaks, we can determine with 
precision the distance between the two merging clumps at each red- 
shift. The merger occurs between z; n = 0.250 and Zf[ n = 0.150. In 
order to follow in detail this event, we pick A^ap = 100 simulation 
snapshots equidistant in redshift between z; n and Zf, n - The redshift 
interval between two consecutive snapshots therefore is Az = 0.001. 
In a ACDM model, this corresponds to a time interval of approxi- 
mately At ~ 10 Myr. 



At redshift z ~ 0.3, the virial mass of our numerical halo is 
~7x 10 A -1 ^/©. In order to increase the lensing efficiency and 
thus to reduce uncertainties in the numerically determined arc cross 
sections, we artificially rescale the cluster mass by multiplying the 
particle masses by a factor / = 2.5. Recalling that the virial radius 
R Y i T of a halo scales as R Y [ T °c where M is the virial mass, 

to obtain a halo of mass / x M which is dynamically stable, we 
also need to rescale the distances by a factor jf 1 ' 3 . This means that 
while the three-dimensional density p remains fixed, the halo sur- 
face density is enhanced by a factor f 1 ^. For this reason we expect 
that, increasing its mass, the numerical cluster will become much 
more able to produce strong lensing effects. 



3.3 Lensing simulations 

Using each of the A snap snapshots, we perform ray-tracing 

simulations . Our method is described in detail elsewhere 

iMeneghetti et alj200dr200lLl2003albT) . but some parameters differ 
from the simulations described earlier. 

For each snapshot, we select again from the simulation box 
a cube of 3 /j^'Mpc side length, containing the high-density re- 
gion of the cluster. The centre of the cube is defined such that the 
selected region contains both merging clumps. A third smaller sub- 
structure enters the field at z ~ 0.23. In the lensing simulations per- 
formed using the snapshots at higher redshift the effect of this small 
clump of matter is not taken into account, because it is out of the 
considered region. We checked its influence on the cluster strong 
lensing properties by selecting a larger region. We found that the 
effect of this substructure is negligible when its distance from the 
main clump is larger than ~ 1.5/7~'Mpc. 

We then determine the three-dimensional density field p of the 
cluster from the particle positions, by interpolating the mass den- 
sity within a re gular grid of 256 3 cells, using the Triangular Shaped 
Cloud method iHocknev & Eastwood 1988). We then produce sur- 
face density fields E by projecting p along three coordinate axes. 
We chose our reference frame such that one axis is perpendicu- 
lar to the direction of merging. By projecting along this coordinate 
axis, we minimize the impact parameter of the infalling substruc- 
ture with respect to the main cluster clump. Hereafter this projec- 
tion will be called "optimal". For comparison, we also investigate 
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Figure 7. Critical lines of the numerically simulated galaxy cluster ("optimal" projection) at the same redshifts shown in Fig.|5] The scale of each panel is 
375". 

the projection along a second axis, where the minimal distance be- 
tween the merging clumps is never smaller than ~ 250 /j^'kpc. The 
third projection, which is not interesting for the purpose of this pa- 
per, since the substructure moves almost exactly along the line of 
sight, will be neglected in the following analysis. 

The lensing simulations are performed by tracing a bundle 
of 2048 x 2048 light rays through the central quarter of the lens 
plane and computing for each of them the deflection angle. Then, a 
large number of sources is distributed on the source plane. We place 
this plane at redshift z s = 1. Keeping all sources at the same red- 
shift is an approximation justified for the purposes of the present 
case study, but the rece nt detect ions of arcs in high- redshift clus- 
ters IZaritskv & Gonzalezl2003tlGladders et alj20 03) indicate that 
more detailed simulations will have to account for a wide source 
redshift distribution. 

The sources are elliptical with axis ratios randomly drawn 
from [0.5, 1]. Their equivalent diameter (the diameter of the circle 
enclosing the same area of the source) is r e = 1". Unlike previous 
studies, we do not distribute the elliptical sources on the central 
quarter of the source plane, since we need to investigate a region 
large enough to contain the caustics of both the merging clumps. 
Nevertheless, we keep the same spatial resolution of the source 
grids in the previous simulations, because sources are initially dis- 
tributed on a regular grid of 64 x 64 instead of 32 x 32 cells. Then, 
we adaptively increase the source number density in the high mag- 
nification regions of the source plane by adding sources on sub- 
grids whose resolution is increased towards the lens caustics. This 
increases the probability of producing long arcs and thus the nu- 
merical efficiency of the method. In order to compensate for this 



artificial source-density increase, we assign for the following sta- 
tistical analysis to each image a statistical weight proportional to 
the area of the grid cell on which the source was placed. We re- 
fine the source grid four times, and the total number of sources is 
typically ~ 30,000. 

Using the ray-tracing technique, we reconstruct the images of 
background galaxies and measure their length and width. Our tech- 
ni que for im age detection and c lassification was described in detail 
bv iBar telmann & Wejsj (l99 4j) and adopted by j Bartelman n et alJ 
( 1998) and lMeneghetti et alJ <200dEoOlL 12003 Jbh . It results in a 
catalogue of simulated images which is subsequently analysed sta- 
tistically. 



4 RESULTS 

In the following sections we show how the lensing properties of our 
numerical model change during the merger. We focus on the "op- 
timal" projection of the numerical cluster, i.e. the projection along 
the axis perpendicular to the direction of the merger, where the ef- 
fects of the merger are expected to be strongest. We refer to the 
second projection only in passing. 



4.1 Critical lines and caustics 

As expected, the critical lines and the caustics of the numerical lens 
evolve strongly during the merger. 

First, we discuss the results for the "optimal" projection. We 
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Figure 8. Maps of the magnification on the source plane produced by the numerically simulated galaxy cluster at the same redshifts shown in Fig.|5] The scale 
of each panel is 1.5 h~ Mpc. 



show the critical lines at some relevant redshifts in Fig.Q At red- 
shift z = 0.230, the main cluster clump and the infalling substruc- 
ture develop separate critical lines. The largest mass concentration 
also produces a small radial critical line (enclosed by the more ex- 
tended tangential critical line). As the merger proceeds, the tan- 
gential critical lines are stretched towards each other. As discussed 
in Sect. 12.21 this is due to the increasing shear in the region be- 
tween the mass concentrations. The critical lines merge approxi- 
mately at redshift z = 0.214. After that, there exists a single critical 
line, which, after a short phase of shrinking, expands isotropically 
while the two clumps overlap. This happens at z ~ 0.203. 

Note that the radial critical line grows during the merging 
event, reaching the maximum extension when the clumps are ex- 
actly aligned and the surface density is highest. Indeed, in order to 
develop a radial critical line, the lens must reach a sufficiently high 
central surface density. 

When the substructure moves to the opposite side of the main 
cluster body, the tangential critical line stretches again and reaches 
its maximum elongation at z ~ 0.190. Then, separate critical lines 
appear around each clump. Their size decreases for decreasing z 
because both the shear and the convergence between the two mass 
concentration decrease as their distance grows. Similarly, the radial 
critical lines shrink. 

During the merger phase, the magnification pattern on the 
source plane changes as shown in Fig. [8] The highest magnifica- 
tions are reached during the phase of maximum overlap, but the 
extent of the highest magnification regions is largest at redshifts 
z = 0.214 and z = 0.190. Therefore, at these redshifts even the 
caustics are most elongated. In Fig. [5] we show the caustics at the 



same redshifts as in Fig. [8] Their evolution reflects that of the crit- 
ical lines. Before redshift z = 0.214, two separated caustics exist; 
as the distance between the two clumps shrinks, their elongation 
in the direction of merging grows. Then the caustics merge into a 
single line, shrink along their long axis and widen perpendicular 
to it. When the two mass concentrations overlap, it has a diamond 
shape with four pronounced cusps. Finally, when the distance be- 
tween the substructure and the main cluster clump grows again, the 
caustic shrinks, elongates in the direction of relative motion of the 
substructure, and finally splits into two small separate caustics. 

In the second projection we have considered, the substructure 
does not pass through the cluster centre. In this case its distance 
from the main cluster clump is always >250/7~'kpc. We show in 
Fig. 1101 the critical lines for this projection. As in the case of the 
"optimal" projection, while the merger proceeds, the critical lines 
of the two merging clumps are stretched towards each other by 
the effect of the increasing shear in the region between the two 
mass concentrations. The largest elongation is reached at redshift 
~ 0.198. However, the substructure does not cross the region en- 
closed by the critical lines of the main cluster clump and there- 
fore the critical lines do not shrink in the direction of merging and 
never expand isotropically. After reaching the maximal elongation, 
it shrinks while the distance between the two clumps grows. The 
caustics (not shown) reflect the same evolution. 



4.2 Tangential and radial arcs 

Since the extent of the high-magnification regions and caustics 
changes during the infall of the substructure onto the main cluster 
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Figure 9. Caustics of the numerically simulated galaxy cluster at the same redshifts as shown in Fig.|5] The scale of each panel is ~ 375" . Note the appearance 
of thin, long caustic structures which are not seen in the analytic model shown in Fig.|3]due to the lack of internal asymmetries and external shear. 



clump as shown in the previous section, we expect that the lens- 
ing efficiency for producing both tangential and radial arcs changes 
accordingly. We quantify the influence of merging on the cluster's 
strong lensing efficiency by measuring its cross sections for differ- 
ent arc properties. 

First, we consider tangential arcs, which can be identified 
among the lensed images due to their large length-to-width ratio 
L/W. By definition, the lensing cross section is the area on the 
source plane where a source must be located in order to be imaged 
as an arc with the specified property. As explained in Sect. 13.31 each 
source is taken to represent a fraction of the source plane. We assign 
a statistical weight of unity to the sources which are placed on the 
sub-grid with the highest resolution. These cells have area A. The 
lensing cross section is then measured by counting the statistical 
weights of the sources whose images satisfy a specified property. If 
a source has multiple images with the required characteristics, its 
statistical weight is multiplied by the number of such images. Thus, 
the formula for computing cross sections for arcs with a property 
p, O p , is 



Willi . 



(19) 



where n, is the number of images of the j-th source satisfying the 
required conditions, and w,- is the statistical weight of the source. 
Using this method, we compute the lensing cross sections for arcs 
with L/W > 5,7.5 and 10, respectively. 

The results are shown in the first three panels of Fig. ll II The 
solid lines refer to the "optimal" projection. All curves exhibit the 
same redshift evolution. Their main properties can be summarized 
as follows: 



• The cross sections grow by a factor of two between z ~ 0.240 
and z ~ 0.220, as shown in Fig.|6| Then, the cross section further 
increases by a factor of five between z ~ 0.220 and z ~ 0.200, i.e. 
within ~ 0.2 Gyr. 

• The curves have three peaks, located at redshifts z\ = 0.214, 
Z2 = 0.203 and Z3 = 0. 190. The peaks at z\ and Z3 correspond to the 
maximum extent of caustics and critical curves along the merging 
direction before and after the moment when the merging clumps 
overlap; the peak at Z2 occurs when the distance of the infalling 
substructure from the merging clump is minimal. 

• Two local minima arise between the three maxima at redshifts 
Z4 = 0.21 1 and zg = 0. 199, where the cross sections are a factor of 
two smaller than at the peaks. At these redshifts, the critical lines 
have shrunk along the merging direction and the caustics are less 
cuspy than at Z2- 

• The cross sections reduce by more than one order of magni- 
tude after z = 0.190. At redshift z ~ 0.180, when the distance be- 
tween the merging clumps is > 1.5 h~ l Mpc, the cross section sizes 
are comparable to those before z ~ 0.240. 

Therefore, during the merger, our simulated cluster becomes 
extremely more efficient in producing tangential arcs. The infalling 
substructure starts affecting the cross sections for long and thin 
arcs when its distance from the main cluster clump is approxi- 
mately equivalent to the cluster virial radius (~ 1.5 h~ 1 Mpc), and 
the largest effects are seen at three different times: 

• When the critical lines (and the corresponding caustics) 
merge, i.e. when the shear between the mass concentrations is suf- 
ficient to produce the largest elongation of the critical lines along 
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Figure 10. Critical lines of the numerically simulated galaxy cluster (second projection) at the same redshifts shown in Fig.|5] The scale of each panel is 375". 



the direction of merging. This happens ~ 100 Myr before and after 
the substructure crosses the cluster centre; 

• when the two clumps overlap, i.e. when the projected surface 
density or convergence is maximal, producing the largest isotropic 
expansion of the critical lines and of the caustics. 

A different behaviour is found in the second projection. The 
results are shown again in the first three panels of Fig. ll li as dashed 
lines. In this case the critical lines and caustics reach a maximal 
elongation at ze = 0.198 and shrink at lower redshifts. Again the 
lensing cross sections reflect the evolution of the critical curves and 
caustics: they reach a maximum extent at id and then their size de- 
creases. A very important result is that, even if the merging clumps 
never get closer than ~ 250 /i~'kpc, the cross sections still grow by 
roughly a factor of five within approximately 50 Myr. 

We now consider the effects of merging on the cross section 
for radial arcs. These are identified from the comple te sample of 
distorted images using the technique described in Meneahetti et al. 
feOOll) . It consists in selecting those arcs for which the measured 
radial magnification at their position exceeds a given threshold. 

The cross section for this type of arcs as a function of red- 
shift is shown in the fourth panel of Fig. 1111 In both projections, 
at redshifts z>0.22, the cross section for radial arcs keeps constant 
and fluctuates around 10 -3 5 /i _2 Mpc 2 . Then, in the "optimal" pro- 
jection, it grows by a factor of five, reaching the highest value at 
z — 0.203. This is the same redshift where the cross sections for 
tangential arcs peak. The enhancement of the convergence due the 
overlapping of the merging clumps thus makes the cluster substan- 
tially more efficient for producing radial arcs. 

Then, the cross section rapidly drops to zero for smaller red- 



shifts. Note that the redshift interval where our cluster is very 
efficient for producing radial arcs is quite small (Az ra( j <0.04). 
Radial arcs have s o far been rep orted in only five galaxy 



Kadtal arcs have so tar been reported in only five galaxy 
clusters ( MS 2137. iFort et alj \l9% : A 370. ISmailetalJll996t 
MS Q440.lGioia et alll998L AC lU. llSfataraian et alll998t A 383, 



ISmith et"all200ll) . Our results contribute to explaining why radial 
arcs are so rare: if they form preferentially during the crossing of 
large substructures through the cluster centre, when the central lens 
surface density is higher, the visibility window of such events is 
very narrow ( < 100 Myr per merger). 

Note that, in the second projection (dashed line), the cross sec- 
tion for radial arcs has some strong fluctuations during the merger 
but never reaches very high values. As for the "optimal" projection, 
it peaks at the same redshift where the cross section for tangential 
arcs is largest, but it changes by less than a factor of two with re- 
spect to redshifts z>0.21. This confirms that only large enhance- 
ments of the central surface density produce a substantial increase 
of the cluster's efficiency for forming radial arcs. 

As mentioned before, assuming a single source redshift is 
justified for the purpose of this work. Indeed, the critical surface 
density changes by ~ 10 per cent when moving the sources from 
redshift one to redshift two, while keeping the lens redshift at 
z\ ~ 0.2. Thus the extent of the critical lines and caustics (and con- 
sequently of the lensing cross sections) is not expected to change 
very much for sources at redshifts above unity. Nevertheless, we 
show in Fig. II li the cross sections obtained by placing the source 
plane at z s = 2 for some characteristic snapshots of our simulations. 
Filled and open circles refer to the "optimal" and the second projec- 
tions, respectively. As expected the relative change in the amplitude 
of the cross sections is modest (less than a factor of 2) for both tan- 
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Figure 11. Lensing cross sections for tangential and radial arcs as functions of redshift. Top left panel: arcs with length-to-with ratio L/W > 5, 05; top right 
panel: arcs with length-to-with ratio L/W > 7.5, O7.5; bottom left panel: arcs with length-to-with ratio L/W > 10, o~io; bottom right panel: radial arcs, a ra( i- 
The solid lines refer to the "optimal" projection, i.e. where the substructure crosses the centre of the main cluster clump, while the dashed lines refer to the 
second projection of the cluster; in both cases the sources are placed at redshift z s = 1. Note the three peaks in the top and bottom left panels; their origin is 
discussed in the text. Cross sections computed considering sources placed at z s = 2 are shown for the "optimal" and second projections by filled and open 
circles, respectively. 



gential and radial arcs. Moreover, the amplitude of the fluctuations 
induced by the merger event on the lensing cross sections seems to 
depend very weakly on the source redshift. 

5 OBSERVATIONAL IMPLICATIONS 

Our results show that cluster mergers could play an important role 
for arc statistics. In particular, since the lensing efficiency grows by 



one order of magnitude during mergers, they may offer a solution 
for the arc statistics problem. 

It is quite important to notice that mergers might have 
some other important observational implications to account for. In 
fact the largest samp le of clusters used for arc statistics studies 
iLuppino et alj[l99^) was selected in the X-ray band, where the 
luminosity is due to bremsstrahlung emission. This is very sensi- 
tive to the dynamical processes going on in the cluster, since it is 
proportional to the square gas density. Therefore, we expect that 





the cluster X-ray luminosity has large variations during a merging 
phase. 

In Fig.HH we show simulated X-ray observations of our clus- 
ter. These were obtained by using X-M AS, a code for simu lating 
data taken with the Chandra satellite l GardinietalJl2003h . The 
cluster emissivity is calculated integrating over the gas density and 
temperature distribution within the cluster simulation, adopting the 



MEKAL plasma model iKaastra & Mewdll993l) with a metallic- 
ity of O.3Z . The upper panels present the results in the "optimal" 
projection at six different redshifts. The upper left panel shows the 
cluster at z = 0.333, before the merger starts. The second panel 
shows the cluster at z = 0.250, when the virial regions merge. The 
third, fourth and fifth panels show the cluster just before, at and just 
after the maximum overlap of the merging clumps, respectively. Fi- 
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Figure 13. X-ray luminosity L x of the numerically simulated cluster as 
function of redshift. The curve is accurately sampled in the redshift range 
z = 0. 15 4- 0.25. The X-ray luminosity has been measured even at z = 0. 1 
and z = 0.333, where is indicated by the filled circles. 

nally, in the sixth panel, the cluster is observed after the end of the 
merger. Since the colour scale is the same for all images, it can be 
easily seen that the X-ray luminosity is highest when the merging 
clumps are closest. Simular results for the second projection are 
shown in the bottom panels. 

In Fig. El we show the observed X-ray luminosity of the nu- 
merical cluster as function of redshift. The curve has a narrow and 
almost symmetric peak located at z ~ 0.200. The X-ray luminos- 
ity grows by more than a factor of four between z ~ 0.300 and 
z ~ 0.200, by roughly a factor of ~ 2.5 between z ~ 0.230 and 
z ~ 0.200 and by roughly a factor of ~ 1.55 between z ~ 0.210 and 
z ~ 0.200. The width at half maximum of the peak is approximately 
Az<0.05. Note that the X-ray luminosity peak is wider than the 
maxima in the arc cross sections, thus the X-ray emission increases 
earlier and decreases later than the arc cross sections during the 
merger. 

If a cluster sample is built by collecting all the objects with 
X-ray luminosity Lx larger than a given threshold, we expect that 
many merging clusters are present among them, since they are 
stronger X-ray emitters. Since these are all very efficient clusters 
for producing gravitational arcs, this could introduce a bias in the 
observationally determined frequency of giant arcs, which could 
become larger than predicted by previous numerical lensing sim- 
ulations in ACDM model. However, it is quite difficult to make 
more robust conclusions here since only one cluster has been anal- 
ysed. Further investigations are needed here. In any case, our re- 
sults show that much caution must be applied when selecting clus- 
ters for arc statistics studies thr ough their X-ray emission (cf. also 
iBartelmann & Steinmetzll996l) . 



6 SUMMARY AND CONCLUSIONS 

In this paper we have investigated how the lensing properties of 
a galaxy cluster change during merging events. Similar dynamical 
processes were not resolved in the previous lensing simulations but 
they might play a relevant role for determining the strong lensing 
efficiency of cluster lenses. 

We address this problem first by using analytic models. When 
simulating a collision between spherical haloes with NFW den- 
sity profiles, we find that both the critical lines and the caustics 



of the lens system strongly evolve during the merger. This be- 
haviour is explained by the change of the shear and of the conver- 
gence induced by the infalling clump. Indeed, while the distance 
between the merging mass concentrations decreases, the shear in- 
tensity grows in the region between the halo centres. The individual 
critical lines and caustics of the main cluster clump and the infalling 
substructure are stretched towards each other until they merge. 

To obtain a quantitative estimate of the impact of mergers on 
the lensing cross sections, we have studied the strong lensing prop- 
erties of a numerically simulated galaxy cluster. Within this numer- 
ical model, a massive substructure falls onto the main cluster clump 
between z; n = 0.250 and Zf, n = 0.150. We have studied the merger 
in detail, picking a large number of simulation snapshots within 
this redshift range. The time separation between consecutive snap- 
shots is approximately ~ 0.01 Gyr. Two different projections of the 
cluster where analysed: in the "optimal" projection, the substruc- 
ture passes through the centre of the main cluster clump, while in 
the second projection the distance between the two mass concen- 
trations is always larger than ~ 250 /z~'kpc. 

The main results of this study can be summarized as follows: 

• As expected from the results of the analytic tests, the shapes 
of critical lines and caustics substantially change during the merger. 
At the beginning, the two clumps develop individual critical lines 
and caustics. These are stretched towards each other while the dis- 
tance between the mass concentrations decreases, because the in- 
tensity of the shear field grows in the region between the approach- 
ing clumps. In the "optimal" projection, when they merge, the re- 
sulting single critical line and caustic shrink along the merging 
direction and then expand isotropically, because of the increasing 
convergence. The same behaviour is observed when the substruc- 
ture moves far away from the main cluster clump after crossing its 
centre. In the other projection the maximum extent of the critical 
lines is reached when the distance between the mass concentrations 
is such that the effect of the shear is largest. After that, the size of 
the critical lines drops. 

• In the "optimal" projection, the lensing cross sections for tan- 
gential arcs change by one order of magnitude during the merger. 
The effect of the infalling substructure starts to be relevant when 
its distance from the main cluster clump is ~ 1.5 /z~'Mpc. The 
cross sections have three peaks located at the redshifts where the 
critical lines have the largest extent along the merging direction, 
or when the shear effect induced by the infalling substructure is 
largest, and at the redshift where the two clumps overlap and conse- 
quently the maximum convergence is reached. Although the effects 
of the merger on the lensing cross sections are important within a 
time interval of ~ 1 Gyr, the strongest impact is thus observed dur- 
ing the central part of the merging phase, on a time scale of ~ 200 
Myr. In the second projection, the lensing cross section for long and 
thin arcs change by a factor of five within a time interval of ~ 100 
Myr. 

• The numerical cluster is highly efficient in producing radial 
arcs only during the merger. The cross section for this type of arcs 
has only one peak, located at the redshift were the infalling sub- 
structure crosses the centre of the main cluster clump. 

Thus, our results show that mergers have a strong impact on 
the strong lensing efficiency of galaxy clusters. Since the lensing 
cross sections for long and thin arcs change by one order of mag- 
nitude during the mergers, these dynamical processes could be a 
possible solution to the arc statistics problem. 

This picture is in principle supported by the fact that samples 
of clusters used in arc statistics studies are selected through their X- 
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ray luminosity, which is very sensitive to the dynamical processes 
arising in the cluster. In particular, we expect that many merging 
clusters are present in these samples, since they are strong X-ray 
emitters. For example, Randall et al. (2002) estimate that the num- 
ber of clusters with luminosities Lx > 5 x 10 44 /i~ 2 erg/sec can be 
increased by a factor of 8.9 due to merger boosts. 

In addition, by surveying cluster s in the LCDCS and in the 
RCS. IZaritskv & Gonzalez! J2003h and lGladders et al] <2003h have 
recently found a high incidence of giant arcs in clusters at high red- 
shift. Their results are particularly interesting since a large number 
of clusters merging at high redshift are predicted by the commonly 
accept ed theory of structure formation. In particular, Gladde rs et alJ 
(2003) speculate that a subset of clusters with low mass and large 
arc cross sections may be responsible for large numbers of arcs in 
distant clusters. One possibility is that such "super-lenses" are clus- 
ters in the process of merging. 

Detailed conclusions are, however, pending on further studies 
quantifying the frequency of mergers and the dependence on arc 
cross sections on the detailed merger parameters, such as the impact 
parameter of the collision, the mass ratio of the merging haloes and 

so forth. Such studies a re now under way. 

In recent studies, fWambsgan ss et all 1200 3l) and iDalal et ail 

(2003) suggest that arc statistics in ACDM models can be recon- 
ciled with the observed abundance of gravitational arcs by adopting 
a broader distribution of source redshifts. Certainly, cross sections 
can change substantially for a cluster of a given mass depending 
on its dynamical state, which makes earlier and current statements 
about the theoretical expectations highly insecure. Moreover, nu- 
merous observational effects need to be taken into account in addi- 
tion for a reliable comparison between numerical simulations and 
observations. 
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